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15.  AC*f'a»l 


A comparison  of  stresses  calculated  using  the  ABS/DAISY 
system  with  those  measured  on  board  the  SL-7  Containership  is 
undertaken  to  verify  the  analytical  procedures  used  in  as- 
sessing the  strength  of  ships  in  a seaway.  The  comparisons 
and  evaluations  are  performed  for  four  different  and  pro- 
gressively more  severe  technical  conditions:  dockside  cali- 
bration, RMS  stresses  in  head  seas,  instantaneous  stresses 
in  head  seas  and  instantaneous  stresses  in  oblique  seas. 

The  overall  comparison  between  calculated  and  measured 
stresses  for  the  dockside  calibration  is  good  where  thermal 
effects  were  small  but  inconclusive  elsewhere.  The  com- 
parison of  RMS  stresses  in  head  seas  is  generally  satisfac- 
tory, using  both  the  spectrum  analysis  approach  and  the 
equivalent  regular  wave  approach,  and  the  comparison  of  in- 
stantaneous stresses  in  head  seas  and  in  oblique  seas  is 
also  good  for  the  wave  conditions  considered. 

The  results  show  that  the  existing  analytical  tools  for 
predicting  wave  loads  and  structural  responses  are  suitable 
to  assess  the  overall  strength  of  the  hull-girder.  All 
the  measured  and  calculated  hull-girder  stresses  are  of  low 
magnitude,  and  no  modifications  to- the  present  hull-girder 
strength  standard  are  deemed  necessary. 
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INTRODUCTION 


This  report  summarizes  the  work  performed  by  the  Re- 
search and  Development  Department  of  the  American  Bureau  of 
Shipping  under  U.S.  Coast  Guard  Contract  No.  DOT-CG-63176-A. 
This  project  is  jointly  sponsored  by  the  Ship  Structure  Com- 
mittee and  the  American  Bureau  of  Shipping,  and  is  part  of 
the  comprehensive  SL-7  Containership  Research  Program  of  ship 
loads,  model  testing,  structural  analysis,  response  analysis, 
full-scale  measurements  and  data  correlation. 

The  long-range  objective  of  the  SL-7  Research  Program 
is  to  advance  the  understanding  of  the  responses  of  hull 
structures  at  sea  and  to  verify  the  current  design  criteria. 
Emphasis  is  focused  on  the  correlation  of  experimental  (full- 
scale  and  model)  data  and  theoretical  predictions  at  different 
phases  of  the  program.  The  overall  research  plan,  historical 
background  and  the  interface  of  various  phases  of  the  program 
were  presented  and  discussed  in  Reference  [1]*.  The  SL-7 
Research  Program  was  sponsored  by  the  Ship  Structure  Committee, 
the  American  Bureau  of  Shipping  and  Sea-Land  Services,  Inc. 

The  major  phases  which  have  been  completed,  either  directly 
under  the  SL-7  program  or  associated  with  it,  can  be  summar- 
ized as  follows: 

1.  Full-scale  Instrumentation 

To  investigate  the  structural  responses  of  a high-speed 
open-deck  containership  at  sea,  an  extensive  full-scale 
instrumentation  system  was  installed  on  board  the  SL-7  class 
containership  S.S.  SEA-LAND  McLEAN  to  measure  wave  heights, 
ship  motions,  accelerations  and  wave-induced  stresses  in  many 
critical  areas.  The  detailed  information  for  this  system  is 
given  in  Reference  [2] . In  addition,  a micro-wave  radar  was 
developed  and  installed  to  measure  wave  elevations.  A finite 
element  analysis  of  the  entire  ship  [3] , using  the  ABS/DAISY 
system  of  computer  programs  [4],  and  a steel  structural  model 
test  [5]  were  carried  out  at  the  planning  and  installation 
stage  of  the  instrumentation  program.  The  results  of  these 
two  research  projects  were  utilized  to  identify  critical 
regions  for  strain- gage  locations. 

After  the  instrumentation  installation  was  completed, 
a dockside  calibration  was  carried  out  by  Teledyne  Materials 
Research  Company  and  reported  in  Reference  [6] . Subsequently, 
a large  amount  of  stress  data  has  been  acquired  for  three 


* Numbers  in  brackets  designate  References  listed  at  the  end 
of  the  report. 
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consecutive  winter  seasons,  between  September  1972  and  March 
1975.  Some  sample  results  are  presented  in  Reference  (7J. 

The  wave-meter  data  was  analyzed  by  Dalzell  (8]. 

2.  Correlation  of  Calculated  Structural  Responses  and  Steel 
Model  Experimental  Data 

Although  the  finite  element  method  utilized  in  predicting 
structural  responses  to  quasi-static  loads  has  been  well  tested 
and  verified  with  full-scale  and  model  experimental  results  in 
recent  years,  it  is  still  desirable  to  validate  the  analysis 
procedure  and  modeling  techniques  in  dealing  with  a special 
structure,  such  as  an  open-deck  containership.  Accordingly, 
the  American  Bureau  of  Shipping  has  performed  a structural 
analysis  of  the  SL-7  Steel  Model,  using  the  ABS/DAISY  system. 

The  calculated  results  together  with  a comparison  with  experi- 
mental data  were  presented  by  Elbatouti,  Jan  and  Stiansen  [9]. 

The  predicted  hull-girder  responses  to  both  bending  and  tor- 
sional loads  were  generally  found  to  be  in  good  agreement  with 
the  measured  data.  Consequently,  the  ABS/DAISY  system  and  the 
current  modeling  practice  employed  at  ABS  are  considered  satis- 
factory in  dealing  with  containerships . 

3.  Correlation  of  Ship-Motion  Calculations  and  Model  Experiments 

Another  significant  phase  of  the  SL-7  program  was  to  verify 
the  wave -load  prediction.  Under  the  sponsorship  of  the  Ship 
Structure  Committee,  a ship-motion  computer  program,  SCORES, 
was  developed  in  1972  based  on  the  two-dimensional  strip  theory 
[10].  It  was  generally  understood  that  the  strip  theory  was 
valid  for  full-form  vessels,  but  its  application  to  fine-form 
ships  was  questionable. 

To  assess  ship  motions  and  wave  loads,  an  SL-7  model  was 
tested  in  oblique  seas  in  1974  [11].  The  first  comparison  of 
the  model  test  data  with  the  SCORES  results  [12]  showed  signi- 
ficant discrepancies  between  the  RAOs  (Response  Amplitude 
Operators)  of  ship  motions  and  wave  loads.  Subsequently,  Oceanics 
Incorporated  introduced  a speed  correction  factor  in  the  existing 
SCORES  program  and  obtained  good  agreement  of  the  RAO  between 
the  model  experimental  data  and  the  theoretical  predictions  [13] . 

Based  on  the  results  of  the  above  three  phases  of  the  SL-7 
research  program,  it  can  be  seen  that  the  finite  element  techni- 
ques and  wave-load  predictions  have  been  generally  validated. 

The  final  phase  in  the  verification  of  the  analytical  procedure 
in  assessing  ship  strength  is  a correlation  of  dynamic  stresses 
of  ships  at  sea. 
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The  objective  of  this  phase  of  the  SL-7  program  is  to 
compare  the  stresses  calculated  using  the  ABS/DAISY  system 
to  those  measured  on  the  SL-7,  in  corresponding  sea  and 
dockside  conditions,  and  to  evaluate  the  results  through  each 
of  four  different  and  progressively  more  severe  technical 
conditions.  Accordingly,  the  work  was  divided  into  the 
following  four  tasks : 

TASK  I - Comparison  of  the  results  of  the  DAISY  stress 
analysis  with  the  results  of  the  full-scale  dockside 
calibration. 

TASK  II  - Comparison  of  the  stress  spectra  calculated 
from  the  DAISY  analysis  results  with  the  selected  full- 
scale  at-sea  stress  spectra. 

TASK  III  - Comparison  of  the  DAISY  stress  analysis  results, 
using  measured  acceleration  and  a specific,  selected  wave 
profile  estimated  from  the  wave  radar  data  in  head  seas, 
with  the  instantaneous  measured  stresses  taken  simultan- 
eously while  the  wave  profile  was  developed. 

TASK  IV  - Investigation  of  TASK  III  for  oblique  sea 
conditions . 
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ANALYSIS  PROCEDURE 


The  procedure  utilized  in  the  performance  of  this  study 
relies  heavily  on  the  ABS/DAISY  system  of  computer  programs. 
A brief  description  of  the  system  and  the  interface  of  the 
associated  element  programs  are  presented  in  Appendix  A. 

The  analysis  comprises  the  following  four  steps: 

1.  Selection  of  Record  Intervals  and  Acquisition  of  the 

Wave  and  Full-Scale  Stress  Data 

To  compare  the  measured  stress  data  with  theoretical 
predictions  for  a ship  at  sea,  it  is  essential  to  have  reli- 
able information  about  the  actual  wave  environment.  Two 
wave-measuring  meters  were  installed  on-board  the  SEA-LAND 
McLEAN,  a Tucker  wave  meter  and  a micro-wave  radar.  A cor- 
relation and  verification  study  [8]  of  the  wave-meter  data 
shows  that  the  Tucker  meter  data  is  in  error  for  the  high 
ship  speeds  of  interest,  and  the  radar  wave  data  appears  to 
be  more  realistic.  Since  there  was  no  meaningful  stress 
data  recorded  with  the  ship  in  a "hove-to"  condition  or  at 
"near  zero"  speed,  the  Tucker  wave  meter  data  could  not  be 
directly  utilized  for  this  study.  Consequently,  the  wave 
environment  was  solely  determined  based  on  recorded  signals 
of  the  wave  radar,  as  presented  by  Dalzell  [8].  However, 
it  should  be  noted  that  the  radar  wave  data  has  not  yet 
been  fully  verified.  Any  possible  errors  incurred  in  the 
wave-measuring  system  or  in  the  data  reduction  procedure 
would  be  essentially  carried  through  the  analysis  and  might 
cause  stress  deviations  in  the  calculation. 

The  selection  of  record  intervals  was  based  on  the 
availability  of  the  measured  stress  data,  the  relative  wave 
angle  and  the  relative  directions  of  the  wave  and  the  swell. 

Regarding  the  full-scale  stress  data,  the  RMS  values, 
stress  spectra  and  stress-time  histories  of  the  midship 
average  bending  sensor  (LVB)  are  presented  in  Reference  [8] 
for  all  the  selected  record  intervals  at  sea.  For  other 
strain  gages,  the  required  stress  data  was  reduced  by  Tele- 
dyne. 


For  Task  I (dockside  calibration) , all  the  test  condi- 
tions were  taken  into  consideration.  The  detailed  information 
about  loading  conditions,  gage  locations  and  the  measured 
data  are  presented  in  Reference  [6] . 
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2.  Approximation  of  the  Wave  Environment  and  Prediction 
of  Wave  Loads 

For  each  selected  record  interval,  the  wave  environ- 
ment was  approximated  either  by  an  equivalent  regular  wave 
or  a wave-spectrum  approach,  depending  upon  the  methods 
utilized  for  each  individual  task.  Details  are  discussed 
under  Tasks  II  and  III. 

Once  the  wave  environment  was  determined,  the  ABS/ 
SHIPMOTION  and  ABS/DYNPRE  computer  programs  were  used  to 
predict  ship  motions,  wave  loads  and  hydrodynamic  pressures 
for  input  to  the  finite  element  structural  model.  The  ABS/ 
SHIPMOTION  program  is  a revised  version  of  SCORES.  ABS/ 

DYNPRE  which  is  an  extended  subroutine  of  the  SHIPMOTION 
program  calculates  the  hydrodynamic  pressure  distribution 
on  the  wet  surface  of  a ship’s  hull. 

As  mentioned  in  the  introduction,  the  revised  version 
of  SCORES,  modified  with  a speed  correction  factor,  has 
been  validated  with  a model  test.  This  modified  version 
is  not  available  to  the  public.  In  the  present  study,  the 
ABS  version  of  the  SCORES  program  (designated  ABS/SHIPMOTION) 
was  utilized  for  wave -load  predictions.  A speed  correction 
factor  similar  to  that  introduced  by  Oceanics  has  been  in- 
corporated into  the  ABS  version.  A sample  comparison  of  the 
RAO  with  the  model  test  data  is  shown  in  Figures  A-2  and  A-3 
of  Appendix  A,  which  shows  that  the  comparison  of  the  RAO 
for  the  vertical  bending  moments  at  midship  is  very  good. 
However,  the  calculated  RAO  for  the  lateral  bending  moment 
at  midship  are  generally  less  than  those  reduced  from  model 
experiments . 

3.  Generation  of  Finite  Element  Structural  Models  and 
Calculation  of  Structural  Responses 

In  calculating  the  hull-girder  responses,  the  entire  ship 
was  first  represented  by  a three-dimensional  coarse-mesh 
finite  element  model.  Subsequently,  fine-mesh  three-dimen- 
sional models  were  utilized  to  determine  the  stress  distribution 
at  selected  strain-gage  locations.  The  ABS/DAISY  computer 
program  was  used  for  this  step  of  the  analysis.  The  struc- 
tural models  and  the  boundary  conditions  used  in  the  analysis 
are  discussed  in  the  following  section. 
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4. 


Analysis  of  Calculated  Results  and  Comparison  with 
Measured  Data 


w 


After  the  completion  of  the  fine-mesh  DAISY  runs  for 
each  task,  a selective  output  containing  strains  and  stresses 
for  the  selected  strain  gage  locations  was  printed  out. 

The  analysis  of  the  calculated  results  varies  task  by  task. 
The  detailed  analysis  methods  and  the  comparison  with  the 
measured  data  are  discussed  separately  under  each  individual 
task. 
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FINITE  ELEMENT  STRUCTURAL  MODELS 


The  SL-7  is  an  880'-6"  x 105'-6"  x 68*-6"  twin-screw 
containership  with  a displacement  of  50,315  long  tons  at 
34'-0"  draft.  A general  arrangement  is  shown  in  Figure  1. 

The  structural  analysis  for  the  SL-7  Containership  was 
performed  using  the  ABS/DAISY  system  of  finite  element  com- 
puter programs  (see  Appendix  A)  for  one  coarse-mesh  model 
comprising  the  entire  ship  and  six  fine-mesh  models  repre- 
senting detailed  structures  of  selected  strain  gage  regions. 


H 
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FIGURE  1 - SL-7  GENERAL  ARRANGEMENT  (FROM  REFERENCE  [3]) 


Three-dimensional  Coarse-Mesh  Model 

The  coarse-mesh  model  consists  of  2602  nodal  points  and 
9866  degrees  of  freedom.  The  major  structural  members  are 
represented  by  rods,  beams,  membrane  plates  and  bending  plates, 


with  a total  of  7122  elements.  The  model  comprised  only  the 
port  side  of  the  vessel  because  of  structural  symmetry  about 
the  centerline  plane.  A three-dimensional  isometric  plot  of 
this  model  is  shown  in  Figure  2.  Computer  plots  of  the  decks, 
bottom,  side  shell  and  the  centerline  profile  are  shown  in 
Figure  3.  Typical  transverse  bulkheads  and  web  frames  in  the 
model  are  shown  in  Figure  4.  For  asymmetric  loadings,  the 
loads  are  divided  into  symmetric  and  anti-symmetric  loading 
conditions,  as  shown  in  Figure  5,  with  appropriate  boundary 
restraints  at  the  centerline  plane,  as  shown  in  Figure  6. 

In  order  to  prevent  rigid  body  movements  due  to  the  possible 
unbalanced  forces  on  the  model,  certain  additional  boundary 
restraints  must  be  imposed  on  the  model,  as  shown  in  Figure  6. 
LOADER,  one  of  the  "DAISY"  preprocessor  programs,  automati- 
cally redistributes  the  unbalanced  forces  in  the  model, 
therefore  minimizing  the  effect  of  these  local  boundary 
restraints  on  the  accuracy  of  the  results. 

Three-dimensional  Fine-Mesh  Models 

Six  fine-mesh  models  in  the  vicinity  of  selected  sensor 
locations,  as  shown  in  Figure  7,  are  utilized  to  determine 
stress  distributions  and  to  facilitate  comparisons  between  the 
calculated  results  and  those  obtained  from  measurements. 

Computer  plots  for  these  six  models  are  shown  in  Figures  8 
through  13,  together  with  sensor  locations  and  designations  as 
specified  in  Reference  16] . Detailed  information  about  all  the 
sensors  installed  on  board  the  S.S.  SEA-LAND  McLEAN  is  given  in 
Appendix  B.  Each  fine-mesh  model  consists  of  rods,  beams  and 
bending  plates,  and  has  the  characteristics  described  in  Table  1. 
Displacements  obtained  from  the  coarse-mesh  analysis  are  used 
as  boundary  conditions  in  the  analysis  of  the  fine-mesh  models. 
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TABLE  1 CHARACTERISTICS  OF  FINE  MESH  FINITE  ELEMENT  MODELS 


Model 

Designation 

Mode  1 

Location 

Number  of 
Nodes 

Number  of 
Elements 

Number  of  Degrees 
of  Freedom 

FK1 

FR. 182-190 
(above  str.no. 1) 

314 

602 

1507 

FM2 

FR. 190-198 

432 

802 

2030 

rM3 

FR. 182-190 
(below  str.no.l) 

391 

536 

1565 

FM4 

FR. 218-238 

318 

680 

1764 

FM5 

FR. 140-150 

326 

522 

1350 

FM6 

FR. 282-298 

393 

541 

1405 

(SENSOR  J6.J7.JI 
( ( AR IA.JB.  JC) 

stnsor  jj.jm.js 

( AN.'A  ,38 , 2C ) ( STBD) 
r-CNSOR  JU.J1.I2 

IAK1A,1B,IC)(  FORT) 


MAIN  DECK 


SECOND  DECK 


TRANE  ISO 
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TASK  I - DOCKSIDE  CALIBRATION 


1.  Loading  Conditions 

The  following  six  calibration  loading  conditions  have 
been  analyzed: 

Condition  1 

Dockside  initial  loading  condition,  with  all  cargo  hold 
and  deck  containers,  except  holds  beneath  Hatches  3,  10  and  14. 

Condition  3 


Deck  containers  removed  from  Hatches  1 through  4 and 
12  through  15. 

Condition  4 

Remaining  deck  containers  on  Hatches  5 through  11  removed. 
Condition  5 


Approximately  one-half  of  containers  removed  from  starboard 
side  of  Hatches  1 through  7 and  from  the  port  side  of  Hatches  8 
through  15,  generating  a torsional  moment.  Hatch  covers  placed 
asymmetrically  to  contribute  to  the  torsional  moment. 

Condition  6 


Completion  of  unloading  described  in  Condition  5.  This 
represents  the  maximum  torsional  load. 

Condition  7 


Nominally  empty  ship  except  for  one  propeller  (47  long 
tons)  loaded  into  Hatch  3 and  one  propeller  in  Hatch  4,  all 
hatch  covers  on. 

Loading  conditions  1,  3,  4 and  7 are  symmetric  about  the 
centerline  plane;  loading  conditions  5 and  6 are  asymmetric 
and  include  torsional  loadings. 


Note;  The  designated  loading  cases  are  identical  with  those 
used  in  Reference  [6] . Loading  condition  2 was  not 
used  during  the  calibration  test. 


17 


The  “SHIPMOM"  program  was  used  to  calculate  static  bending 
moments  and  to  generate  hydrostatic  pressures  for  the  DAISY 
model.  These  SHIPMOM  results  were  compared  with  the  results  of 
a previous  less  refined  SHIPMOM  computer  run  with  fewer  stations, 
described  in  Reference  [6).  The  comparison  indicated  very  little 
difference  between  the  two  runs. 

2.  Comparison  of  Calculated  Stresses  and  Measured  Data  - 
Longitudinal  Stresses 

The  comparison  of  stresses  measured  by  strain  gages  and 
those  calculated  at  corresponding  gage  locations  is  presented  in 
Tables  2 and  3.  In  Table  2 the  comparison  was  made  by  subtrac- 
ting loading  condition  1 from  each  loading  condition,  taking 
loading  condition  1 as  a datum  loading.  Similarly,  Table  3 
takes  loading  condition  4 as  the  datum  loading. 

A preliminary  investigation  of  the  experimental  and 
analytical  stress  results  indicates  some  instances  of  agreement, 
and  some  instances  of  disagreement.  The  correlation  of  full- 
scale  measurements  and  computed  results  can  only  be  established 
after  carefully  verifying  the  data  and  the  environmental  condi- 
tions . 

The  experimental  data  for  midship  sensors  1 ( LVB ) , 15  (LSTS) 
and  18  (LSTP) , where  LVB  measures  the  average  longitudinal  verti- 
cal bending  stress  of  port  and  starboard,  LSTP  and  LSTS  measure 
the  longitudinal  stress  at  top,  port  and  starboard  respectively, 
show  some  discrepancies. 

Sensor  number  1 should  average  the  values  of  the  longitu- 
dinal stress  components  of  the  LSTS  and  LSTP  sensors.  However, 
differences  in  the  comparative  stress  values  appear  with  the 
change  in  the  reference  of  datum  loading  condition,  or  the 
so-called  zero-stress  reference  case.  Table  4 illustrates 
these  differences  when  referring  loading  condition  7 to 
loading  conditions  1,  3 or  4 . 

Loading  condition  1 represents  no  change  in  mechanical 
loading  of  the  ship  from  the  initial  calibration  condition  (zero 
reading  for  the  strain  gages)  that  took  place  when  the  ship  was 
navigating  through  the  Maas  river,  but  the  measured  stresses 
range  between  -2677  and  1654  psi  (Reference  [6)). 

Loading  condition  7 is  predominantly  a symmetric  loading 
case,  but  the  measured  stresses  from  the  symmetrically  located 
sensors  15  and  18  indicate  asymmetric  response,  especially  in 
loading  condition  7-4. 
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TABLE  2 COMPARISON  OF  CALCULATED  AND  MEASURED  STRESSES  (PSI) 
USING  LOADING  CONDITION  1 AS  DATUM  LOADING 


; 96 (TGMS3) 
97 (TGMS4 ) 
98 (TGMS1X) 
99 (TGMS2X) 
ICO (TGMS3X) 
101 (TGMS4X) 
(17<LSBS) 


[20(LSBP) 
73(R11A) 
76 (R12A) 
79 (R13A) 
82 (R14A) 
30  (ARIA) 
33 (AR2A) 


36 (AR3A) 
39  (AR4A) 


43 (R1A) 
46 (R2A) 

4 

49 (R3A) 
52 (R4A) 


170  376 

-123  -1882 
531  1785 

531  589 

-85  -115 

-133  142 

-525  -1096 
-454  -1864 


-719  -593  -2 


-719  152  -2 

-1362  -1084  -4 
-1212  -1964  -4 
85  -937  1 

85  -725  1 

98  -513  2 


111  -163 


1803  -365 

863  1150 

-267  1099 

-802  589 

-10  399 

515  687 

2420  3257 

1300  -195 

1476  1491 

-1297  -1125 
-1884  -1939 
266  1145 

-1136  -2418 
1974  948 

1386  343 

-246  -1254 
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TABLE  3 COMPARISON  OF  CALCULATED  AND  MEASURED  STRESSES  (PSI) 
USING  LOADING  CONDITION  4 AS  DATUM  LOADING 


Loading 
Cond . 


Sensor 


(5-4) 

(6-4) 

(7-4) 

Calc. 

Meas. 

Calc. 

Meas. 

Calc. 

Meas , 

786 

707 

1249 

1016 

2643 

2252 

687 

-991 

1000 

1306 

2643 

3738 

885 

870 

1497 

1374 

2643 

1877 

-1377 

-2042 

-2059 

-3034 

-3 

117 

-1287 

-3419 

-1876 

-2427 

183 

-1986 

1098 

342 

1720 

798 

89 

-799 

1812 

3214 

2803 

4229 

761 

958 

-1316 

-2101 

-1958 

-2568 

-416 

-642 

-27 

-192 

33 

-1103 

486 

-55 

1079 

-457 

1694 

0 

-133 

-628 

375 

789 

617 

902 

738 

846 

-251 

-687 

-283 

-549 

-1162 

-1465 

-469 

-181 

-818 

136 

-1162 

-227 

-43 

228 

-112 

742 

384 

1141 

324 

764 

437 

764 

593 

873 

1531 

2577 

2079 

3181 

2524 

3071 

760 

-1339 

999 

837 

2089 

1896 

1310 

1426 

1753 

1700 

784 

-549 

-698 

-2184 

-1020 

-2717 

784 

-1926 

-953 

-2395 

-1412 

-3821 

1475 

-970 

623 

2345 

718 

2127 

1350 

1582 

-738 

-1536 

-1306 

-2084 

540 

-768 

980 

1952 

1803 

2510 

540 

1898 

633 

1825 

1184 

2168 

612 

2024 

-270 

-763 

-462 

-872 

672 

I-  'll 

20 
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It  can  be  seen  that  the  degree  of  agreement  between  the 
measured  and  calculated  results  varies  significantly  from 
Table  2 to  Table  3.  This  indicates  some  inconsistencies  of 
the  measured  data  between  different  loading  conditions,  probably 
due  to  changes  in  ambient  temperatures. 

According  to  Table  IV  of  Reference  [6] , a temperature 
difference  of  15°  F between  port  and  starboard  sides  was  recor- 
ded during  the  calibration  test.  A temperature  difference  of 
21°  F between  the  deck  and  sea  water  was  also  recorded.  With 
these  magnitudes  of  temperature  gradient,  the  thermal  stress 
could  be  as  high  as  1500  psi.  However,  due  to  the  lack  of 
sufficient  information  about  the  temperature  distribution 
during  the  calibration  test,  it  is  impossible  to  incorporate 
thermal  stresses  into  the  analysis. 

To  evaluate  the  effects  of  these  temperature  changes, 
longitudinal  stresses  measured  at  17  selected  sensors  are 
shown  in  Figures  14  through  16.  The  sensors  are  described  in 
Table  5. 

Figure  14  shows  the  longitudinal  stresses  for  calibration 
loading  conditions  3 through  7,  utilizing  loading  condition  1 
as  the  datum  loading.  The  deviations  of  the  measured  stresses 
from  the  calculated  values  are  generally  within  a bandwidth  of 
+ 1,500  psi.  Based  on  the  temperature  differential  recorded 
(Turing  the  calibration  test,  shown  in  Table  6,  the  maximum 
thermal  stress  may  be  as  high  as  1,500  psi.  With  this  thermal 
stress  margin  in  mind,  the  overall  comparison  of  the  measured 
and  computed  stresses  shown  in  Figure  14  is  reasonably  good. 

In  an  attempt  to  minimize  the  possible  thermal  effects, 
changes  in  stresses  between  two  loading  conditions  with  the 
least  temperature  differential  were  also  examined.  The  results 
of  loading  case  4-3,  which  represents  pure  vertical  bending, 
are  shown  in  Figure  15,  With  a few  exceptions,  the  stress 
deviations  fall  within  a bandwidth  of  + 400  psi.  This  magni- 
tude of  deviation  is  regarded  as  acceptable,  considering  the 
sensibility  and  reliability  of  strain-gage  readings.  The 
results  of  loading  case  7-6,  which  reflects  both  vertical 
bending  and  torsion,  are  shown  in  Figure  16.  With  the  exception 
of  Sensors  30  and  39,  the  stress  comparison  is  generally  good. 

To  minimize  the  possible  thermal  effects,  it  is  advisable, 
for  the  future  calibration  tests,  that  other  means  of  loading 
be  used  to  create  appreciable  mechanical  strains  and  that 
complete  steel  temperature  data  of  the  deck  and  side  shell 
(port  and  starboard)  be  recorded. 
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TABLE  4 MEASUREMENTS  OF  LONGITUDINAL  STRESSES  (PSI)  IN 
MIDSHIP  SENSORS 


Sensor  No. 


Loading  Condition 


1S(LSTS) 

18 ( LSTP) 

Average  of  15  and  18 
l(LVB) 


(7-1) 

(7-3) 

(7-4) 

2252 

4008 

3738 

1789 

3158 

1877 

2020.5 

3583 

2807.5 

2119 

3267 

2252 

TABLE  5 SELECTED  SENSOR  LIST 


SENSOR 

NUMBER 


SENSOR 
NOMEN . 


SIGNAL  NOMENCLATURE 


Longitudinal  Vertical  Bending 
Longitudinal  Stress  Top  Starboard 
Longitudinal  Stress  Top  Port 
Longitudinal  Stress  Bottom  Starboard 
Longitudinal  Stress  Bottom  Fort 

R • Foreward  Rosettes  on  the 
Main  Deck,  See  Figure  11 


AR  *=  Ai't  Rosettes  on  the  Main 
Deck,  See  Figure  1? 


R «=  Forward  Rosettes  on  the 
Main  Deck,  See  Figure  13 


Another  uncertainty  in  comparing  the  measured  and  calcu- 
lated stresses  is  the  influence  of  plate  unfairness  and  local 
plate  bendings.  So  far  as  the  longitudinal  stresses  in  the 
hull  structures  are  concerned,  this  kind  of  influence  is 
considered  insignificant.  All  the  longitudinal  strain  gages 
were  installed  on  very  thick  plates  where  the  effect  of  unfair- 
ness, if  any,  would  not  be  generally  noticeable.  Where  local 
plate  bending  was  apparent,  such  as  the  bottom  hull  plating, 
the  strain  gages  were  located  at  the  quarter  span  between 
frames  (floors)  to  minimize  local  plate  bending.  In  addition, 
all  the  primary  plating  was  represented  by  bending  plate 
elements  in  the  fine-mesh  models.  Panel  bendings  are  also 
included  in  the  calculation.  Consequently,  the  plate  unfairness 
and  local  plate  bending  should  not  be  regarded  as  an  influential 
parameter  to  the  stress  discrepancies  in  this  case. 


Sensor  Number 


Calculated 
Stress  (jjp  1 ) 


♦ 3 0 0 0 


Notes : 

Mots  arc  based  or.  results  of 
Lo.ulir.v  Condition.*  3-1, 4-1,  S- 
6-1  a:.d  7-1.  ac  listed  in 
Table  2 


FIGURE  14  - DEVIATION  OF  MEASURED  LONGITUDINAL  STRESSES  FROM  CALCULATED  VALUES 


Sensor  Humber 


'leasmed 
Stress  (psi) 


♦2000 


Calculated 


Calculated 
Stress  (pal) 


r 


1 


3.  Comparison  of  Calculated  Stresses  and  Measured  Data  - 
Transverse  Stresses 

The  calculated  and  measured  transverse  stresses  in  the 
transverse  box  girder  at  Fr.  194-196  are  plotted  in  Figure  17 
for  loading  cases  3-1,  4-1,  5-1,  6-1  and  7-1.  Regarding  the 
difference  between  the  measured  and  calculated  stresses  at 
each  gage  location,  these  figures  show  approximately  the  same 
results  as  those  shown  in  Table  2.  However,  by  comparing  the 
stress  distribution  patterns,  it  can  be  seen  that  the  measured 
data  reflects  generally  a higher  degree  of  bending  and  torsion 
loads,  probably  caused  by  thermal  expansions. 

To  minimize  the  thermal  effect,  loading  cases  4-3  and  7-6 
were  again  selected  for  comparison.  The  results  are  shown  in 
Figure  18.  No  improvement  in  stress  agreement  can  be  gained 
for  loading  case  4-3,  which  represents  a pure  longitudinal 
bending  load. 


LOADING  CASE  3-1 


LOADING  CASE  1-1 


LOAPING  CASE  >-I 


COMPRESSION 


\ 


INBOARD 


o MEASURED 


-FORWARD 


SCALE 


TENSION 


0 X 2 
KM  I 


FIGURE  1?  - CAPARISON  OF  CALCULATED  AND  MEASURED  TRANSVERSE  STRISSES  IE  THE 

TRANSVERSE  BOX-GaROER  AT  I*R.  19H-1SI 
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Since  the  primary  loads  of  the  transverse  box  girder  are 
induced  by  hull-girder  bending  and  twisting,  in  addition  to 
thermal  loads,  manufacturing  imperfections  could  also  contri- 
bute to  the  deviations  of  the  measured  data.  Furthermore, 
the  strain  gages  have  a single  element  which  measures  the 
strain  in  the  transverse  direction  only.  The  measured  stress 
data  does  not  account  for  the  Poisson's  effect  which  may  cause 
significant  stress  variations.  However,  a reasonably  good 
agreement,  particularly  of  the  stress  distribution  patterns, 
can  be  seen  for  loading  case  7-6,  with  the  only  exception  of 
Sensor  95. 

In  conclusion,  a reasonably  good  agreement  of  the  measured 
stresses  and  the  calculated  values  has  been  obtained  for  several 
selected  calibration  conditions  by  minimizing  the  possible 
thermal  effects. 
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For  Task  II,  a stress  comparison  between  the  calculated 
and  measured  RMS  values  in  head  seas  was  carried  out  by  an 
equivalent  regular  wave  approach  and  by  a spectrum  analysis, 
based  on  a study  of  the  relationship  between  computed  wave- 
induced  stresses  and  wave  heights. 


In  order  to  examine  this  relationship,  a regular  wave 
with  a length  of  808.5  feet  and  with  its  crest  at  the  midship 
was  selected  for  ship  motion  and  stress  calculations.  Three 
different  wave  heights  as  shown  in  Table  7 were  taken  into 
consideration,  with  a shir  speed  of  10.9  knots  in  head  seas. 


TABLE  7 WAVE  CHARACTERISTICS  AND  SHIP  MOTION  DATA  - LINEARITY  STUDY 


Notes: 

a All  conditions  consist  of  full  cargo  loads  and  head  waves 
b Heave  is  positive  down 
c Pitch  is  positive  bow  up 


The  ship's  motions  and  accelerations  which  are  required  as 
part  of  the  finite  element  model  input  were  calculated  using  the 
SHIPMOTION  computer  program.  A set  of  quasi-static  pressures 
(still-water  plus  wave  profile)  and  inertia  forces  corresponding 
to  the  calculated  accelerations  were  applied  on  the  finite  ele- 
ment model  and  a correction  was  made  to  ensure  that  the  applied 
vertical  shearing  forces  and  bending  moments  were  generally 
identical  to  those  obtained  directly  from  the  SHIPMOTION  program. 
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The  calculated  results  are  shown  in  Figures  19,  20,  and  21 
at  three  selected  locations  along  the  length  of  the  ship.  As 
shown  in  Figure  19,  the  wave-induced  longitudinal  and  transverse 
stresses  in  the  main  deck  plating  at  the  midship  section  are 
approximately  in  direct  proportion  to  the  wave  height.  They 
depart  from  linearity  at  stations  remote  from  the  midship  section 
(Figures  20  and  21) , even  though  the  hull- girder  shearing  forces 
and  bending  moments  vary  linearly  with  wave  heights,  as  predic- 
ted by  a ship  motion  calculation. 

The  non-linearity  is  due  to  the  combination  of  the  local 
bending  in  the  transverse  direction  and  the  wave-induced  forces 
in  the  longitudinal  direction. 

Recent  model  experiments  in  towing  tanks  indicate  that  the 
motions  and  wave  loads  do  not  vary  linearly  with  wave  heights. 

A realistic  trend  cannot  yet  be  defined  due  to  the  lack  of 
sufficient  experimental  data.  Based  on  the  results  shown  in 
Figures  20  and  21,  it  can  be  seen  that  the  hull- girder  responses 
(wave-induced  stresses)  exhibit  a non-linear  pattern  even  with 
linear  shear  and  bending  loads. 


WAVE  1IEICHT.  METERS 

riGURE  19  - WAVE-INDUCED  STRESSES  VS  WAVE  HEIGHTS  - 

mat;.’  deck  tlatijig  at  the  midship  section 


In  the  equivalent  regular  wave  method,  a given  wave  spec- 
trum can  be  approximated  by  a regular  wave  with  equivalent 
energy.  With  this  approach,  no  assumption  concerning  the 
relation  between  stresses  and  wave  heights  is  necessary.  On 
the  other  hand,  this  equivalent  regular  wave  approach  will  not 
account  for  the  variation  in  responses  to  different  wave  fre- 
quencies . 

A spectrum  analysis,  which  depends  completely  upon  an 
assumption  of  linearity  between  stresses  and  wave  heights, 
is  valid  for  performing  stress  comparisons  only  in  the  midship 
region. 

After  a thorough  review  of  these  two  methods,  it  was 
decided  to  perform  the  RMS  stress  comparison  using  both 
approaches  selectively.  Three  recorded  intervals  were  selec- 
ted for  the  equivalent  regular  wave  approach.  For  comparison, 
the  same  intervals  were  also  used  for  the  spectrum  analysis. 
Because  of  the  non-linear  effects,  the  spectrum  analysis  was 
carried  out  for  the  longitudinal  stresses  at  midship  only. 

The  stress  RAO  (Response  Amplitude  Operator)  generated  for 
these  three  intervals  was  then  utilized  to  calculate  the 
stress  spectrum  for  other  selected  intervals. 

The  detailed  procedures  of  these  two  approaches  and  the 
computed  results  are  discussed  in  the  following  sections. 

1.  Method  1 - Equivalent  Regular  Wave  Approach 

A.  Description  of  Method 


a.  Calculation  of  Equivalent  Regular  Waves 

The  equivalent  wave  encounter  frequency,  wave  height 
and  wave  length  of  the  regular  waves  were  calculated  based  on 
the  radar  wave  spectrum  shown  in  Reference  [8] , using  the 
following  formulas 


r“e2 

j^ei  du>£ 

Equivalent  wave  encounter  frequency  = 


(1) 


( UJe2 

Uel  V“e>  d“e 


Equivalent  wave  height  = 2.5 


e2 


el 


W d“c 


1/2 


(2) 
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Equivalent  wave  length  = 


(3) 


•1 


2irg 

~^2 


u is  derived  from  the  equation 
u cosB  o 

ID  = (D y=F  (D* 


where  g = gravity  constant 


u * Wave  frequency 

o)e  * encounter  frequency,  measured  directly 
from  the  time  history 


8 = heading  angle 


(4) 


V«e) 

Wel'Ue2 


ship's  speed 

wave  spectrum  in  the  frequency  domain 

lower  and  upper  bounds  of  frequency 
within  the  domain  under  consideration 


Three  wave  conditions  were  selected  for  this  task.  The 
detailed  information  and  the  equivalent  waves  for  these  three 
wave  conditions  are  shown  in  Table  8. 


TABLE  8 WAVE  CONDITIONS  SELECTED  FROM  REFERENCE  (8]  FOR 
COMPARISON  OF  RMS  STRESSES  IN  HEAD  SEAS 


wave 

Condition 

Tape 

Index 

Interval 

Run  No. 

Ship  Speed 
(knots) 

Equivalent 
Wave  Length 
(ft.) 

Equivalent 
Wave  Height 
(ft.) 

4 

145 

18 

5 

405 

10.9 

808.5 

20.92 

5 

145 

24 

29 

429 

18.7 

808.5 

21.97 

6 

145 

29 

50 

450 

28.1 

561.5 

16.47 
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b.  Calculation  of  Wave-Induced  Loads 


A computer  subroutine  was  developed  by  ABS  for  the  calcu- 
lation of  hydrodynamic  pressures  (ABS/DYNPRE)  acting  on  the 
ship  in  a seaway.  However,  it  was  found  that  the  wave-induced 
vertical  bending  moments  obtained  from  an  integration  of  the 
hydrodynamic  pressures  were  not  in  full  agreement  with  those 
obtained  from  the  ABS/SHIPMOTION  program.  In  order  to  elim- 
inate such  differences,  the  following  two  modifications  were 
utilized  to  correct  the  vertical  bending  moments  obtained  by 
integration  of  the  imposed  pressures  and  inertial  forces. 

Method  1A  - Quasi-static  Pressures 

A quasi-static  pressure  approach  was  used  to  calculate 
the  input  load  for  the  ABS/DAISY  runs.  In  this  approach,  the 
quasi-static  pressures  and  the  inertia  forces  at  the  instan- 
taneous ship  position  as  determined  by  the  ship  motion 
calculation  were  first  applied  to  the  structural  model.  Then, 
the  bottom  pressures  were  modified  to  make  the  vertical  wave 
bending  moments  comparable  with  those  obtained  from  the  ship 
motion  program.  This  was  the  procedure  used  in  the  determin- 
ation of  the  relationship  between  computed  wave-induced  stresses 
and  wave  heights  described  at  the  beginning  of  Task  II. 

Method  IB  - Hydrodynamic  Pressures 

The  hydrodynamic  pressure  coefficients  obtained  from  the 
ABS  External  Pressure  subroutine  (DYNPRE)  were  modified  at 
each  ship  station  to  make  the  vertical  wave  bending  moments 
comparable  with  those  obtained  from  the  ship  motion  calculation. 

In  both  the  ship  motion  and  the  external  pressure  calcu- 
lations, information  is  obtained  only  for  regions  of  the  ship 
under  the  still-water  line.  To  account  for  the  actual  wave 
profile,  the  external  pressures  were  linearly  extended  up  or 
deleted  down  to  the  wave  surface,  depending  on  whether  the 
wave  surface  was  above  or  below  the  still-water  line.  By  so 
doing,  the  total  external  pressure,  which  includes  both  the 
hydrodynamic  and  the  static  components,  becomes  zero  at  the 
wave  surface.  The  typical  pressure  distributions  for  wave 
surfaces  above  and  below  the  still-water  line  are  shown  in 
Figures  22  and  23  respectively.  The  wave  characteristics  and 
ship  motion  data  used  for  Methods  1A  and  IB  are  shown  in  Table  9. 

c . Structural  Analysis 

After  the  calculation  of  dynamic  and  static  loads,  the 
stress  analyses  for  coarse-mesh  and  subsequently  for  fine- 
mesh  models  were  performed  by  using  the  ABS/DAISY  finite 
element  program.  For  the  comparison  of  wave-induced  stresses, 
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TABLE  9 WAVE  CHARACTERISTICS  AND  SHIP  MOTION  DATA  - METHOD  X - TASK  II 


Method 

Wave 

Cond. 

(a) 

Load- 

ing 

Cord. 

Ship 

Speed 

(knots) 

Wave  Length 
(ft.) 

Wave  Height 
(Peak-to-Trough) 
(ft.) 

(b) 

Heave 

(ft.) 

(c) 

Pitch 

(deg.) 

Location  of 
Wave  Crest 
Forv.urd  From  A.P. 
(ft.) 

■ 

11 

10.9 

808.5 

20.92 

-0.69 

-2  JO 

220.0 

1A 

12 

10.9 

808.5 

20.92 

3.93 

256 

440.0 

■ 

13 

10.9 

808.5 

20.92 

-0.4l 

1.73 

660.0 

■ 

11 

10.9 

808.5 

20.92 

-OJ59 

-2  JO 

220.0 

IB 

12 

10.9 

808.5 

20.92 

353 

256 

440.0 

■ 

13 

10.9 

808.5 

20.92 

...... 

-OAl 

1.73 

- - . 

660.0 

— 

14 

18.7 

808.5 

21.97 

152 

-1.93 

220.0 

IB 

5 

15 

18.7 

808.5 

21.97 

7.30 

3-27 

440.0 

16 

18.7 

808.5 

21.97 

-2.65 

1.2 

660.0 

17 

28.1 

mm 

16.47 

-a  63 

0.13 

220.0 

IB 

6 

18 

28.1 

1 

16.47 

0.84 

0.01 

440.0 

19 

28.1 

D 

16.47 

-a  68 

-0.15 

660.0 

Notes : 


A,  b,  and  c see  sign  convention  under  Table  7 


the  net  dynamic  responses  can  be  obtained  simply  by  deducting  the 
still-water  response  from  the  overall  results  for  each  wave 
loading  condition. 

d.  Generation  of  Stress-time  History  Curve 

In  order  to  determine  the  amplitude  of  the  wave-induced 
stresses  for  each  wave  condition,  a stress-time  history  of  one 
complete  stress  cycle,  as  shown  in  Figure  24,  was  plotted  based 
on  three  different  wave- crest  positions  at  the  midship  and  the 
quarter  points,  and  each  wave-crest  position  was  treated  as  a 
separate  loading  condition. 
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DAISY  Output 


-IdralireJ  Sties* 

Cut  vc 


Note:  0l,2  or  3 • Strata  Calculated 
b«*»d  on  wavo  crost  located  at 
L/4  from  A.P.,  Midship,  and 
L/4  from  F.P. 


FIGURE  24  - XDEALIIF.D  STRESS-TIME  HISTORY  CURVE,  BASED  ON 
THREE  WAVE  CREST  POSITIONS 


Two  assumptions  were  made  in  generating  the  stress-time 
history.  First,  the  time  variation  of  stress  was  assumed  to 
be  sinusoidal.  Second,  the  calculated  longitudinal  wave- 
induced  stress  was  assumed  to  be  in  phase  with  the  wave-induced 
vertical  bending  moment  at  the  ship  station  under  consideration. 
In  the  comparison  of  the  calculated  and  the  measured  RMS  stres- 
ses, the  double  stress  amplitude  (peak-to-trough)  was  utilized. 

B.  Calculated  Results 

For  comparison,  the  calculated  and  measured  RMS  stresses 
for  the  three  selected  wave  conditions  are  shown  in  Tables  10 
and  11  and  are  also  plotted  in  Figures  25,  26  and  27.  The 
measured  values  were  based  on  data  reduction  performed  by 
Teledyna.  The  RMS  stresses  calculated  by  Methods  1A  and  IB 
are  generally  in  good  agreement,  as  shown  in  Table  10.  These 
stresses  are  generally  of  the  same  order  of  magnitude  as  the 
measured  values.  The  calculated  and  measured  RMS  stresses  for 
sensors  ARIA,  AR2A  and  AR3A  in  Table  11  show  good  agreement; 
however,  sensor  AR4A  consistently  shows  that  the  measured  data 
are  much  lower  than  the  calculated  values. 
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TAoiiE  10  COMPARISON  OF  CALCULATED  AND  MEASURED  RMS  LONGITUDINAL 
STRESSES  (PEAK-TO-TROUGH,  PSI) 


Sensor 

Calc.  RMS 

Meas . RMS 

No, 

Location 

Stress 
(Method  IB) 

Stress 

1 (LVB) 

MN.DK. MID- 

9,673 

6,743 

SHIP  SEC.' 

(AVG.  OF 

STBD.  AND 

PORT) 

8,469 

6,344 

6,887 

5,368 

15(LSTS) 

MN.DK. MID- 
SHIP SEC. 

6,887 

5,840 

(STBD) 

17 (LSBS) 

NEAR- BOTTOM 

PL.  (STBD) 

4,570 

4,164 

3,039 

3,297 

18<LSTP) 

MN.DK. MID- 
SHIP SEC. 

9,673 

5,367 

(PORT) 

6,887 

5,529 

20 (LSBP) 

NEAR- BOTTOM 

4,570 

2,974 

PL.  (PORT) 

4,623 

2,743 

3,039 

2,876 

30 (ARIA) 

MN.DK. FR. 
143-144 

5,285 

5,027 

(PORT) 

5,850 

4,769 

3,433 

4,649 

33 (AR3A) 

MN.DK.FR. 

143-144 

5,285 

5,587 

(STBD) 

5,850 

5,767 

3,433 

4,939 

36 (AR3A) 

MN.DK.  FR. 
143-144 

10,505 

7,537 

(STBD) 

11,657 

7,563 

6,844 

6,576 

39 (AR4A) 

MN.DK.FR. 

143-144 

10,499 

4,333 

(STBD) 

11,450 

4,017 

6,930 

3,786 

2 . Method  2 - Stress  Spectrum  Approach 

A.  Description  of  Method 

The  procedures  used  to  calculate  wave-induced  dynamic 
loads  and  to  perform  the  structural  analysis  are  the  same  as 
those  described  in  the  previous  sections.  For  the  external 
pressures,  Method  IB  was  used  in  conjunction  with  wave  spectra 
from  Reference  [ 8 ] for  the  selected  wave  conditions.  Since  the 
selected  three  basic  intervals  have  three  different  ship  speeds, 
ranging  from  10.9  to  28.1  knots,  a stress  RAO  curve  has  to  be 
generated  for  each  ship  speed.  In  order  to  have  sufficient 
points  for  each  selected  wave  spectrum  a total  of  23  wave 
lengths,  as  shown  in  Table  12,  was  used.  The  wave  lengths 
were  determined  in  such  a way  that  these  three  basic  wave 
spectra  could  be  closely  represented. 

For  each  wave  length  the  relative  wave- crest  position 
was  selected  to  maximize  the  vertical  bending  moment  at  the 
ship  station  under  consideration.  This  information  was  readily 
available  from  the  ship  motion  calculations.  Assuming  that  the 
longitudinal  stresses  in  the  deck  plating  are  in  phase  with 
the  vertical  bending  moment,  the  stress  amplitude  could  be 
directly  calculated  for  each  wave  frequency. 

In  an  attempt  to  minimize  the  possible  stress  deviations 
due  to  the  non-linear  relation  between  stresses  and  wave  heights, 
an  average  wave  height  of  19.68  ft  was  used  for  all  the  differ- 
ent waves.  The  computed  stresses  were  then  divided  by  the  wave 
amplitudes  to  generate  the  stress  RAO.  Once  the  RAO  curve  in 
terms  of  a unit  wave  amplitude  for  the  selected  stress  sensor 
is  constructed,  a stress  spectrum  can  be  obtained  by  multiply- 
ing the  ordinate  of  the  wave  spectrum  by  ( RAO ) ^ . The  measured 
wave  spectrum  has  a low-frequency  cut-off  value,  as  illustrated 
in  Figures  28  through  41,  below  which  the  data  are  unreliable. 


B. 


Calculated  Results 


I 


The  calculated  stress  spectra  for  the  three  selected  basic 
wave  conditions  are  shown  in  Figures  28,  29  and  30  for  sensor 
LVB  and  Figures  31,  32  and  33  for  sensors  LSBS  and  LSBP.  The 
measured  stress  spectra  reproduced  from  Reference  [8]  are  also 
shown  for  comparison.  The  calculated  RMS  stresses,  expressed 
in  terms  of  2/3a  (where  A is  the  area  under  the  stress  spectrum 
curve)  are  listed  in  Table  14. 

In  order  to  make  a further  comparison  for  sensor  LVB,  eight 
additional  recorded  intervals  with  nearly  head  sea  conditions 
and  comparable  ship  speeds  were  selected  from  Reference  [8]. 

The  particulars  for  these  intervals  are  listed  in  Table  13. 
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TABLE  12 


WAVE  CHARACTERISTICS  AND  SHIP  MOTION  DATA  - METHOD  2 - TASK  II 
WAVE  HEIGHT  (PEAK-TO-TROUGH)-  19.68  FT.  (6M) 


*ve 

lond. 


la) 

Load- 

ing 

Cond. 

Ship 

Speed 

(Knots) 

Wave 

Encounter 

Frequency 

Irad./sec. ) 

Wave  Length 
(ft.) 

m 

(c) 

Pitch 
(deg. ) 

Location  of 
Wave  Crest 
Forward  From 
(ft.) 

20 

10.9 

.279 

3366.9 

9.32 

0.19 

471.7 

21 

10.9 

.451 

1460.5 

7.26 

0.81 

448.8 

22 

10.9 

.560 

1015.85 

5.21 

1.42 

439.0 

23 

10.9 

.643 

808.5 

3.69 

1.85 

433-8 

24 

10.9 

.729 

682.96 

2.76 

1.95 

432.0 

25 

10.9 

.816 

550.4 

1.75 

0.96 

436.8 

26 

10.9 

1.088 

350.8 

-0.39 

0.31 

564.3 

27 

18.7 

.430 

1901.8 

8.31 

0.87 

496.8 

28 

18.7 

.548 

1295.4 

7.38 

1.55 

464.3 

29 

18.7 

.611 

1093-14 

6.97 

2. 00 

452.0 

30 

18.7 

.745 

821.5 

6.58 

2.89 

436.7 

31 

18.7 

.866 

644.46 

2.14 

1.97 

431.4 

32 

18.7 

1.019 

509.33 

1.75 

0.02 

431.8 

33 

18.7 

1.138 

434.5 

1.36 

-0.05 

426.8 

34 

28.1 

.305 

3887.9 

7.57 

0.65 

820.3 

35 

28.1 

.508 

1758.7 

8.23 

1.51 

541.2 

36 

28.1 

.479 

1322.0 

8.72 

2.25 

492.4 

37 

28.1 

.712 

1088.0 

9.63 

3.00 

465.6 

38 

28.1 

.817 

880.9 

8.27 

3.81 

438.9 

39 

28.1 

.967 

695.7 

-1.52 

1.30 

433.4 

40 

28.1 

1.130 

561.5 

1.00 

0.05 

427.6 

41 

28.1 

1.273 

47S.4 

1.03 

-0.01 

420.7 

42 

28.1 

1.579 

359.3 

-.14 

0.0 

578.2 

Note* » 

a,  b,  and  c see  sign  convention  under  Table  7 


Since  these  intervals  have  ship  speeds  comparable  with 
wave  condition  6,  it  is  reasonable  to  use  the  stress  RAO  curve 
for  wave  condition  6,  (Figure  30b)  to  calculate  the  stress 
spectra  for  the  additional  intervals.  The  calculated  stress 
spectra  together  with  the  measured  stress  and  wave  spectra 
are  shown  in  Figures  34  through  41.  The  calculated  RMS  stresses 
for  these  additional  intervals,  together  with  the  measured  values, 
are  also  listed  in  Table  14. 


As  shown  in  Figures  28  through  41,  the  calculated  stress 
spectra  generally  exhibit  shapes  similar  to  those  obtained 
from  the  measured  data.  Some  of  the  calculated  peak  frequen- 
cies and  peak  amplitudes  deviate  significantly  from  the 
measured  curves.  However,  the  calculated  RMS  stresses  as 
shown  in  Table  14  generally  agree  well  with  the  measured  values. 


For  a further  comparison  of  the  three  different  methods  to 
calculate  external  pressures,  some  calculated  and  measured  RMS 
longitudinal  stresses  in  the  main  deck  at  midship  (Sensor  LVB) 
are  shown  in  Table  15. 


TABLE  15  COMPARISON  ON  CALCULATED  AND  MEASURED  RMS  VERTICAL  BENDING 
STRESSES  (PEAK-TO-TROUGU,  PSI)  EOR  SENSOR  LVB  AT  MIDSHIP 


B.  As  shown  in  Tables  14  and  15,  the  RMS  stresses  calculated 
by  the  spectrum  analysis  (Method  2)  agree  very  well  with  those 
obtained  from  the  measured  data.  On  the  other  hand,  at  some 
strain-gage  locations  (Table  10)  the  RMS  stresses  computed  by 
Method  1 are  significantly  greater  than  the  measured  values. 

This  may  be  attributed  to  the  fact  that  the  equivalent  regu- 
lar wave  approach  ignores  the  possible  variations  in  structural 
responses  to  different  wave  frequencies.  The  spectrum  analysis, 
which  accounts  for  all  the  significant  wave  frequencies,  gives 

a better  average  value.  It  should  be  noted  that  the  spectrum 
analysis  is  valid  only  where  the  relation  between  the  stresses 
and  wave  lengths  can  be  approximated  as  linear. 

C.  The  calculated  stress  spectra,  as  shown  in  Figures  28 
through  41,  assume  similar  shapes  as  those  obtained  from  the 
measured  data.  However,  the  calculated  peak  frequency  and  peak 
amplitude  may  deviate  from  the  measured  ones  in  some  comparisons. 

D.  The  stresses  calculated  by  Method  1A  are,  as  expected, 
comparable  with  those  obtained  by  Method  IB  for  head  sea 
conditions . 


Text  continued  on  Page  53 
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TASK  III  - INSTANTANEOUS  STRESS  COMPARISON  IN  HEAD  SEAS 

1.  Selection  of  Record  Intervals 

Two  record  intervals  were  selected  for  this  phase  of  the 
study.  The  detailed  environmental  information  on  these  two 
selected  intervals  is  shown  in  Figures  42  and  43,  as  obtained 
from  Reference  18] . The  selected  time  span  for  stress  correl- 
ation covers  two  complete  encounter  cycles,  shown  shaded  within 
a circle  on  the  time  history  of  the  radar  wave  elevation. 
Figures  42  and  43. 

2.  Ship  Motion  Calculation 

To  facilitate  the  wave -load  calculations,  the  recorded 
wave  elevation  within  the  time  period  of  two  complete  cycles 
was  approximated  by  a regular  wave  with  a mean  encounter 
frequency,  an  average  height  and  an  equivalent  wave  length 
(See  the  enlarged  scale  in  Figures  42  and  43) . 

The  wave  particulars  for  the  two  selected  conditions 
are  shown  in  Table  16.  Once  the  wave  particulars  were  esta- 
blished, a ship  motion  calculation  was  performed  using  the 
ABS/SHIPMOTION  program  to  compute  motions,  accelerations, 
shearing  forces  and  bending  moments  for  three  positions  of  the 
wave  crest  along  the  length  of  the  vessel,  namely,  at  the  mid- 
ship and  two  quarter-length  stations. 


TABLE  16  WAVE  CHARACTERISTICS  AND  SHIP  MOTION  DATA  - TASK  III 


Wave 

Cond. 

Load- 

ing 

Cond. 

Ship 

Speed 

(knots) 

Wave  Length 

(ft.) 

Wave  Height 
( Peak-t  o-Trougn ) 

Heave 

(ft.) 

Pitch 
(deg. ) 

43 

10.9 

1551.0 

51.2 

-15.9 

-2.8 

15 

1)4 

10.9 

1651.0 

51.2 

20.4 

1.48 

45 

10.9 

1651.0 

51.2 

11.4 

4.79 

46 

28.1 

808.5 

30.6 

14.35 

0.48 

16 

47 

28.1 

808.5 

30.6 

4.64 

5.40 

48 

28.1 

808.5 

30.6 

-15.64 

-1.99 

Location  of 
Wave  Crest 
Forward  From  A.P. 
(ft.) 


FEET 


DATE  VO  TIME 
POSITION 
COURSE  AND  SrE  :0 
SEA  STA‘C 
HAVE  WEIGHT 
' REE  DIR 
SHELL  HEIGHT 
* REL  DIR 

VISUAL  hEAThCR  f CONSENTS 

OCAST  / 


Cl  - 1 C- 74  2*1  CO 

A3* 29  N 24-51  w 

2-30  . 10.9  KNOTS 


wave  KIOHT  g -ATIST1CS  CEEET) 

tuc«.er/dtn.  kao^r-vdar 

P-T  SAMPLE  SIZE  90  " 70  121 

MAXIMUM  HEIGHT  IS. 9 43  4 S7.2 

I 0Th  HIGHEST  HTS  15.0  39.0  39.7 

3RD  HIGHEST  HTS  13.0  33. D 32. S 


4.0  RMSt  SpCCTRAi  1 13  3 34.0  33.S 


RADAR  HAVE  ELEVATION 


TIME  MINUTES 


FIGURE  42  - CONDITION  15  - VOYAGE  32W  — TAPE  145 
RUN  405 


INDEX  18  --  INTERVAL  5 


LOO  ROOK  DATA 


01- 12-74  1630 

41-07  N 40-09  14 

2- ’S  . 29.1  KNOTS 


DATE  AND  TIME 
POST  TICN 
COURSE  AND  SPEED 
SEa  state 
WAVE  HEIGHT  2S  t’EET 

• REL  DIR  2 STBD 

SHELL  HEIGH-  20  ‘‘EET 

’ REL  DIR  2 STRD 

VI SLAV  HDA'hER  / COMMENTS 


. STATISTICS  (fEETl 
TlCKER/DYM.  HE  AD/RAO  AR 


MAXIMUM  HEIGHT 
10-H  HIGHEST  HTS 
3RD  HIGHEST  HTS 


RADAR  HAVE  ELEVATION 


time  MINUTES 


FIGURE  43  - CONDITION  16  - VOYAGE  32W  — TAPE  145  --  INDEX  29 
RUN  450 
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3.  Input  to  DAISY  Model 

A modification  of  hydrodynamic  pressure  coefficients, 
described  as  Method  IB  in  Task  II,  was  carried  out  to  make 
the  wave-induced  vertical  shearing  forces  and  bending  moments 
acting  on  the  DAISY  model  generally  identical  to  those  calcu- 
lated directly  by  the  SHIPMOTION  program.  A sample  comparison 
is  shown  in  Figure  44. 


;g  i i ■*  i..;*'  ,<■  C'.i^  ..  i.»  . i \c-*  j 

1 1 * • « *•»  C r,t"  'A  t.o  * I ii  l J*  i*  l.v,  '’n* " 1*7*  '•*>  M ' './o  i 


u 7 - 

O } — 

u a — 

S2  -i  — 

- -2  -A . P ■ 
« - 
< -3  — 

^-4U 


-Shearing  Force  — — 


i 


. . • 1 ‘ . . 


i 


oM  iL-.l!  • i*  —i4,  —i,?  1 ~ ..j ••  -i« * _jr' \ ri “ . j i!Z  ‘ c J —j-t-  iV~  1-1^. 

[.  ; |j  FT  m Ir  • . * u>  To  »*«  In  T*vTvAiu  tv*  Ti/  fo  T*4  V-6 


FIGURE  44  - WAVE-INDUCED  BENDING  MOMENTS  AND  SHEARING  FORCES 


4.  Comparison  of  Calculated  Stresses  and  Measured  Data 


For  each  wave  condition,  three  values  of  stress  were 
calculated,  corresponding  to  the  three  wave- crest  positions 
described  above.  In  order  to  construct  a diagram  of  stress- 
time history  based  on  the  three  calculated  values,  it  was 
assumed  that  the  computed  stress  follows  a sinusoidal  curve 
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with  a frequency  equal  to  the  wave -encounter  frequency. 

It  was  further  assumed  that  the  stress  variation  is  directly 
proportional  to  the  dominant  wave -load  component,  the  phase 
angle  of  which  can  be  obtained  from  the  ship  motion  calcula- 
tions. For  longitudinal  stresses  in  the  deck  or  bottom 
regions,  the  dominant  wave -load  component  was  assumed  to  be 
the  vertical  bending  moment.  A mean  stress  amplitude  obtained 
from  these  three  calculated  stress  values,  with  proper  phase 
angles,  was  then  utilized  to  plot  an  idealized  stress  curve 
(See  Figures  45  and  46) . 

A.  Mean  Midship  Vertical  Bending  Stress  - Sensor  IVB 

The  calculated  and  measured  mean  midship  vertical  bending 
stresses  (Sensor  LVB)  are  shown  in  Figures  45  and  46  for  Wave 
Conditions  15  and  16  respectively.  The  measured  stress-time 
history  was  selected  within  the  same  time  period  used  for  the 
selection  of  wave  data.  Recognizing  the  difference  between 
the  recorded  wave  elevation  and  the  idealized  wave  profile, 
it  can  be  seen  that  the  comparison  of  the  calculated  stresses 
and  the  measured  data  is  very  good. 

B.  Other  Stress  Gages  at  Midship  for  Wave  Condition  15 

In  addition  to  the  mean  midship  vertical  bending  stress 
(Sensor  LVB) , the  following  stresses  at  the  midship  were  also 
recorded  at  the  same  instant  as  the  selected  wave  data: 

LSTP  - Longitudinal  Stress  Top  Port 

LSBP  - Longitudinal  Stress  Bottom  Port 

LSBS  - Longitudinal  Stress  Bottom  Starboard 

The  measured  and  calculated  longitudinal  stresses  at  the 
locations  of  Sensors  LSTP  and  LSBP,  as  shown  in  Figures  47  and 
48,  exhibit  a good  agreement.  However,  the  comparison  of 
measured  and  calculated  stresses  at  the  location  of  Sensor  LSBS, 
Figure  49,  is  relatively  poor.  The  measured  stresses  of  LSBP 
and  LSBS  should  be  generally  comparable  in  head  seas.  The 
difference  between  the  measured  stresses  of  LSBP  and  LSBS 
may  be  attributed  to  irregularity  of  the  encounter  wave  and 
.ossibly  irregularity  of  gage  signals. 

It  should  be  noted  that  all  the  stress  data  other  than 
the  mean  midship  vertical  bending  stress  (LVB)  was  generally 
recorded  on  Recorder  No.  2,  whereas  the  wave  data  and  LVB 
were  recorded  on  Recorder  No.  1.  In  addition,  all  the  stress 
data  other  than  LVB  was  reduced  separately  from  that  of  the 
wave  data. 


56 


CP3I 


C.  Three-arm  Rosette  Gages  at  Frame  143-144  for 
Wave  Condition  16 

Other  sensors  recording  simultaneously  with  the  selected 
wave  data  for  Wave  Condition  16  are  rosette  gages  on  the  deck 
at  Frame  143-144.  The  longitudinal  elements  of  these  rosettes 
are  designated  by  ARIA,  AR2A , AR3A  and  AR4A . 

Since  an  element  of  the  rosette  measures  strain  in  one 
direction  only,  the  measured  stress  reflects  a straight  product 
of  the  strain  component  and  Young's  modulus  (E) , but  not  the 
true  stress.  For  purposes  of  comparison,  it  is  essential  to 
calculate  stresses  on  a consistent  basis.  Therefore,  the 
computed  stresses  shown  in  Figures  50-53  are  the  product  of 
the  longitudinal  strain  component  and  Young's  modulus. 

The  comparison  of  the  measured  and  calculated  stresses 
at  locations  of  Sensors  ARIA,  AR2A , AR3A  and  AR4A  as  shown 
in  Figures  50  through  53  are  generally  good. 

5.  Conclusions 

Based  on  the  results  discussed  above,  the  following 
conclusions  can  be  established  for  Task  III: 

A.  The  comparison  of  the  measured  and  calculated  longitudinal 
stresses  for  a selected  period  of  time  in  head  seas  is  gener- 
ally good.  The  calculated  stress  amplitudes  agree  well  with 
the  mean  values  of  the  selected  two  stress  cycles,  with  the 
exception  of  two  gage  locations. 

B.  For  the  selected  cases,  the  wave  data  are  generally 
consistent  with  the  stress  data. 

C.  To  calculate  structural  responses  in  head  seas,  irregular 
waves  can  be  approximated  by  regular  waves  with  equivalent 
amplitude  and  frequency. 
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TIME  MINUTES 


COMPARISON  cr  THE  CALCULATED  AND  MEASURED  INSTANTANEOUS 
LONGITUDINAL  STRESSES  (SENSOR  LSTP) , WAVE  CONDITION  15 


TIME  MINUTES 


COMPARISON  OF  THE  CALCULATED  AND  MEASURED  INSTANTANEOUS 
LONGITUDINAL  STPESSES  (SENSOR  LSBP) , WAVE  CONDITION  15 


FIGURE  48 
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TASK  IV  - INSTANTANEOUS  STRESS  COMPARISON  IN  OBLIQUE  SEAS 


1.  Method  of  Approach 

Due  to  the  complexity  of  the  wave  load  pattern  in  oblique 
seas,  difficulties  were  experienced  in  dealing  with  the  hydro- 
dynamic  pressures  and  balancing  of  the  structural  model.  To 
improve  the  consistency  between  the  external  pressures  and 
the  ship  motion  calculation,  hydrodynamic  pressure  coefficients 
were  modified  for  lateral  moments  in  addition  to  the  modifi- 
cations for  vertical  moments  described  in  Task  II. 

To  investigate  the  feasibility  of  performing  an  instan- 
taneous stress  comparison  in  oblique  seas,  two  record  intervals 
were  selected  for  test  runs.  The  detailed  information  for 
the  selected  intervals  and  the  radar  wave  data  as  obtained 
from  Reference  [8]  are  shown  in  Table  17. 


TABLE  17  WAVE  CONDITION  SELECTED  FROM  REFERENCE  (8]  FOR  THE 

COMPARISON  OF  INSTANTANEOUS  STRESSES  IN  OBLIQUE  SEAS 


l Wave 

I Condition 

Tape 

Index 

Interval 

Run  N'o. 

i 

| 

1 

Ship  Speed 
(knots) 

Wave 

Relative 

Dir. 

(deg. ) 

Heading  Angle 

(deg . ) 

17 

143 

11 

44 

345 

32.3 

64  Port 

116 

18 

143 

12 

48 

349 

31.8 

41  Port 

139 

.,.  — 1 

For  instantaneous  stress  comparison,  two  complete  cycles 
of  the  recorded  wave  elevation  were  selected  and  approximated 
by  a regular  wave  with  a mean  encounter  frequency,  an  average 
height  and  an  equivalent  wave  length,  (see  the  enlarged  scale 
in  Figures  54  and  55) . 

For  each  wave,  three  wave-crest  positions,  at  midship  and 
the  two  quarter-length  points,  were  utilized  to  calculate 
wave  loads.  Each  wave-crest  position  was  treated  as  a separ- 
ate loading  condition.  The  particulars  of  these  loading 
conditions  are  shown  in  Table  18. 
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loo  ao~n  data 


DATE  AND  T'.n:  01  2Ci 

PCSITJCN  | •1S-00  N 13 
CCLf-oL  A\J  SP'.ry  | ?l,7  . 30.3 

SC*  STATE  1 B 
MAVC  nLlOi.T  j 13  fCCT 
RCL  DIR  o*1  PCPT 
SWELL  hCIOnr  10  FEET 
1?-L  O'  I 'I  6A  PCST 

— vtsc-v.  i.:a’h*:b  / co'-.'tMs 


'Measured  Wave 


TlA'iF-W>N.  HEAD,  WAS 


TIME  MINUTES 


FIGURE  54 


CONDITION  1?  - VOYAGE  3 IK  - TAFE  143  - INDEX  U - INTERVAL  44 
RUN  345 


01-C5.7A  2.100 

-15-C9  N 13- 30  « 
266  . 31  .R  WOTS 


-Measured  Wave 
Regular  Wave  . 


15  PCET 
.11  PC'!  I 
1C  FEET 


Notes ; 


a,  b,  and  c see  sign  convention  under  Table  7 

^ Heading  angle  is  measured  counterclockwise  from  ship 
centerline  to  wave  direction 

e Roll  is  positive  starboard  deck  edge  down 


2.  Calculation  of  Wave-induced  Loads 

To  determine  the  wave-induced  loads  for  input  to  DAISY,  the 
SHIPMOTION  and  DYNPRE  programs  were  used.  In  oblique  seas,  the 
wave-induced  moments  consist  generally  of  three  components,  namely 
vertical,  lateral  and  torsional  moments.  Since  the  wave-induced 
loads  obtained  by  integration  of  a set  of  hydrodynamic  pressures 
and  the  corresponding  inertia  forces  would  generally  not  agree 
with  those  obtained  from  the  ship  motion  calculation,  a correction 
coefficient  for  the  hydrodynamic  pressure  was  introduced  to  modify 
the  vertical  and  lateral  moments  for  the  external  pressures.  After 
this  modification,  the  torsional  moments  may  not  be  in  agreement 
with  those  obtained  from  the  ship  motion  calculation. 

Because  of  the  complexity  of  the  wave -load  pattern  in  oblique 
seas  and  the  difficulty  in  balancing  the  DAISY  model  in  the  later- 
al direction,  the  external  pressures  were  not  extended  up  or 
deleted  down  to  the  wave  surface  for  this  task. 
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3.  Structural  Analysis 


The  procedures  utilized  for  calculating  structural  responses 
for  this  task  are  generally  the  same  as  those  in  Tasks  II  and  III, 
except  that  the  asymmetric  wave  load  was  divided  into  symmetric 
and  anti-symmetric  components  which  were  then  treated  as  separate 
loading  conditions  in  DAISY.  The  combination  of  the  symmetric 
and  anti-symmetric  cases  gives  the  total  response  for  both  the 
port  and  starboard  sides.  (See  Figure  5) . 

For  each  wave,  in  addition  to  the  three  dynamic  loading 
conditions  corresponding  to  three  wave-crest  positions,  a static 
case  reflecting  the  still-water  condition  was  also  included  in 
the  analysis.  The  wave-induced  stresses  were  then  obtained  by 
subtracting  the  still-water  stresses  from  the  total  responses 
for  each  dynamic  loading  condition. 

4.  Calculated  Results 


In  order  to  generate  a stress-time  history  curve  based  on 
the  calculated  stresses  at  three  different  wave-crest  positions, 
the  following  two  assumptions  were  made: 

a.  The  curve  was  assumed  to  be  sinusoidal. 


b.  The  calculated  longitudinal  wave-induced  stress  was 
assumed  to  be  in  phase  with  the  effective  moment  (M  ) combining 
the  wave-induced  vertical  and  lateral  bending  moments,  expressed 
in  terms  of 


M 


e 


\ + Ok 

(SMjVsry 


(5) 


where  M , M = wave-induced  vertical  and  lateral  bending 

L moments  respectively,  at  the  location  under 

consideration. 

SMy , SM  = hull-girder  section  modulus  with  respect  to 
vertical  and  lateral  bending  moments  respec- 
tively, at  the  location  under  consideration. 


The  phase  angles  of  and  M were  obtained  directly  from  the 
ship  motion  calculation,  and^the  effective  moment,  M , was  then 
determined  graphically.  A sample  plot  is  shown  in  Figure  56. 
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Bending 


7x10 


FIGURE  56  - GRAPHIC  SOLUTION  FOR  THE  PHASE  ANGLE  OF  THE  EFFECTIVE  BENDING  MOMENT 


5.  Comparison  of  Calculated  Stresses  and  Measured  Data 


For  each  wave  condition,  three  values  of  stress  were  calcu- 
lated, corresponding  to  the  three  wave -crest  positions  described 
above.  In  order  to  construct  a diagram  of  stress-time  history 
based  on  the  three  calculated  values,  it  was  assumed  that  the 
computed  stress  follows  a sinusoidal  curve  with  a frequency 
equal  to  the  wave- encounter  frequency.  It  was  further  assumed 
that  the  stress  variation  is  in  direct  proportion  with  the 
dominant  wave -load  component,  and  that  the  phase  angle  for 
combined  wave-induced  vertical  and  lateral  bending  moments  can 
be  obtained  using  the  method  described  in  the  previous  section. 

A mean  stress  amplitude  obtained  from  these  three  calculated 
stress  values  with  proper  phase  angles  was  then  utilized  to 
plot  idealized  stress  curves  as  shown  in  Figures  57  through  67. 


Based  on  the  results  shown  in  Figures  57  through  67,  it 
can  be  seen  that  the  calculated  stress-time  histories  generally 
agree  with  the  measured  data,  except  for  Sensor  AR4A  in  Figure  67. 
Text  contained  on  page  74. 
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CALCULATED  TOTAL  STRESSES 


In  the  structural  analysis,  a total  of  46 
conditions,  with  wave  heights  up  to  51.2  feet, 
into  consideration.  The  detailed  characterist 
wave  conditions  are  shown  in  Tables  7,  9,  12, 
To  examine  the  overall  structural  responses  to 
wave  conditions,  the  fine-mesh  results  were  sc 
stresses  greater  than  11,360  psi  (800  kg/cm^) . 
selective  stress  outputs  which  exhibit  the  max 
for  all  the  wave -loading  conditions  are  shown 
through  C5  in  Appendix  C.  In  those  tables  all 
shear,  principal  and  Hencky-von  Mises  stresses 
The  printed  values  represent  the  total  stresse 
water  plus  wave)  in  kg/cm^. 
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The  maximum  total  stresses  in  the  selected  fine-mesh 
models  which  repres  nt  the  most  critical  regions  of  the  hull 
girder  are  general  j less  than  20,000  psi  (1408  kg/cm-). 

The  only  exceptio'  a are  the  plate  element  Nos.  329,  333  and 
347  of  fine-mesh  model  5 (Table  C-4  and  Figure  12) , where 
the  highest  stress  reaches  21,513  psi  (1515  kg/cm-).  Most 
of  the  high  stresses  occur  at  the  hatch  corner  region  where 
the  structural  discontinuity  causes  significant  stress  con- 
centration . 


The  strength  deck  of  the  SL-7  class  containership  was 
constructed  of  ABS  H33  higher-strength  steel.  In  accordance 
with  the  current  ABS  Rules,  the  permissible  hull-girder 
bending  stress  for  this  vessel  is  30,100  psi.  Therefore, 
the  highest  calculated  total  stress,  which  includes  stress 
concentration,  is  approximately  71%  of  the  permissible 
hull-girder  bending  stress. 

Based  on  these  calculated  stresses  and  the  full-scale 
stress  data  measured  to  date,  no  stresses  greater  than  the 
permissible  limit  have  been  found. 
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CONCLUSIONS 


1.  The  overall  comparison  between  calculated  and  measured 
stresses  for  the  dockside  calibration  is  generally  inconclusive 
because  of  significant  temperature  differentials  during  the 
test  and  the  low  magnitudes  of  the  applied  loads.  However, 
good  agreement  between  the  calculated  and  measured  stresses 
was  obtained  when  thermal  effects  were  small.  For  future  cali- 
bration tests,  it  is  recommended  that  complete  temperature  data 
be  recorded  and  that  appreciable  mechanical  strains  be  generated 
in  the  structure. 

2.  The  comparison  of  calculated  and  measured  RMS  stresses  in 
head  seas  is  generally  satisfactory,  using  both  the  spectrum 
analysis  approach  and  the  equivalent  regular  wave  approach.  The 
correlation  using  the  former  approach  which  takes  into  account 
the  variation  in  response  to  different  wave  frequencies,  shows 
better  agreement  than  that  obtained  from  the  latter  approach. 
However,  it  should  be  noted  that  the  spectrum  analysis  approach 
is  valid  only  near  the  midship  section  where  the  wave-induced 
stresses  are  generally  in  direct  proportion  to  wave  heights. 

3.  The  agreement  of  calculated  and  measured  instantaneous 
stresses  in  head  seas  is  generally  good;  the  calculated  stress 
amplitudes  agree  well  with  the  mean  values  of  the  measured 
stresses  over  a time-span  of  two  complete  encounter  cycles. 

4.  The  agreement  of  calculated  and  measured  instantaneous 
stresses  in  oblique  seas  is  also  generally  good  for  the  wave 
conditions  considered. 


5.  Based  on  the  results  obtained  from  this  project,  it  can  be 
concluded  that  the  existing  analytical  tools  for  predicting  wave 
loads  and  structural  responses  are  suitable  to  assess  the  overall 
strength  of  the  hull  girder.  All  the  measured  stress  data  reduced 
to  date  from  the  SL-7  instrumentation  program  and  all  the  calcula- 
ted hull-girder  stresses  from  the  present  study  were  found  to  be 
of  low  magnitude.  Consequently,  no  modifications  to  the  present 
hull  girder  strength  standard  are  deemed  necessary. 

6.  Regarding  the  structural  responses,  the  ABS/DAISY  system  in 
its  present  form  is  considered  satisfactory  for  performing  static 
analyses  for  either  quasi-static  or  dynamic  loads. 

7.  In  order  to  assess  the  strength  of  local  structures  and  also 
to  improve  the  current  methodology , further  research  is  deemed 
necessary  to  improve  the  calculation  method  for  external  pressures 
(possibly  verified  by  model  experiments)  and  to  improve  the  ship 
motion  program  (SCORES)  to  account  for  three-dimensional  effects 
and  possibly  for  non-linear  effects  with  respect  to  wave  heights. 
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APPENDIX  A - ABS/DAISY  COMPUTER  PROGRAM  SYSTEM 


The  primary  components  of  the  DAISY  computer  program 
system  are  the  preprocessor  programs,  the  DAISY  finite  element 
program  itself  and  the  postprocessor  programs.  The  flow  chart 
of  the  system  is  shown  in  Figure  A-l. 

Preprocessor  System 

Seven  principal  computer  programs  form  the  nucleus  of 
the  three-dimensional  ship  structure  preprocessing  system. 

Each  program  performs  a specific  function  and  is  used  in  a 
particular  sequence  so  that  the  output  of  one  program  can  be 
used  as  input  to  another.  The  seven  programs  are  listed  below 
in  the  order  of  their  use : 

1.  SHIPMOM  (SHIP  MOMENT)  calculates  hull- girder  shear 
forces,  bending  moments  and  vertical  deflections 
for  a vessel  in  still  water  or  statically  poised  on 
a wave. 

With  the  vessel's  hull  geometry  and  a description 
of  lightship  weight  and  cargo  weights  as  inputs,  the 
equilibrium  draft  and  trim  are  calculated  for  a ship 
in  a still-water  condition  or  in  any  sinusoidal  or 
trochoidal  wave.  The  sinusoidal  waves  are  directly 
applicable  to  ship  dynamics  and  their  effect  on  the 
static  response  of  the  vessel  differs  very  little 
from  that  of  the  trochoidal  waves  traditionally  used. 
The  program  calculates  quasi-static  lateral  bending 
moments  and  torsional  moments  if  the  vessel  is  poised 
in  an  oblique  wave.  Computer  line  plots  of  the  shear 
force  diagrams,  bending  moment  diagrams  and  deflec- 
tions are  automatically  generated. 

2.  SHIPMOTION  (ABS  version  of  SCORES  program)  based  upon 
the  two-dimensional  ship  theory  for  six  degrees  of 
freedom  predicts  a vessel's  motions  (velocity  and 
acceleration)  in  a seaway,  longitudinal  and  lateral 
wave-induced  bending  moments  and  shear  forces,  as 
well  as  torsional  moments  due  to  waves.  The  seaway 
can  be  regular  sinusoidal  waves,  irregular  long-crested 
waves,  or  regular  and  irregular  short-crested  seas. 

The  sea  condition  can  also  be  represented  by  measured 
wave  data.  The  program  computes  the  ship  motions  and 
the  dynamic  components  of  bending  moment  and  shear 
force  resulting  from  the  sea  state  chosen,  as  well  as 
long-term  values,  using  probability  theory  based  on 
the  intended  service  of  the  vessel.  A sample  compari- 
son of  the  RAO  for  the  vertical  and  lateral  bending 
moments  with  model  experimental  data  are  shown  in 
Figures  A-2  and  A-3. 


DYNPRE  (DYNAMIC  PRESSURE)  converts  the  data  gener- 
ated  by  SHIPMOTION  into  dynamic  load  inputs  to  the 
DAISY  finite  element  program.  DYNPRE  converts  the 
hydrodynamic  pressure  distribution  generated  by 
SHIPMOTION  into  a pressure  distribution  at  the 
nodal  points  of  the  finite  element  structural  model 
and  modifies  the  pressures  to  compensate  for  certain 
discrepancies  in  the  SHIPMOTION  hydrodynamic  pres- 
sure calculation. 

EXAM  generates  the  ship  structural  finite  element 
model  for  either  a portion  of  or  for  the  entire  hull 
structure.  Basic  data  generated  by  the  program 
consists  of  nodal  points  and  loadings  applied  to  the 
elements  and  nodes.  The  program  can  generate  a finite 
element  model  of  any  type  of  vessel. 

With  the  vessel's  hull  geometry  finite  element  nodal 
point  locations  are  generated  with  a minimum  of  user 
input.  Structural  properties  of  bending  or  membrane 
plates , beams  and/or  rod  members  are  input  in  a 
convenient  tabular  form  from  which  the  program  auto- 
matically generates  appropriate  elements  to  represent 
the  ship  structure.  The  elements  are  generated  and 
connected  in  a way  that  minimizes  the  bandwidth  of 
the  structural  stiffness  matrix. 

EXAM  automatically  calculates  the  symmetric  and  anti- 
symmetric components  of  a general  seaway  loading. 

Loads  in  the  form  of  pressure  heads  are  automatically 
calculated  and  applied  to  the  elements  according  to 
the  vessel's  hydrostatic  cargo  and  sea  loadings 
described  in  SHIPMOM.  The  user  also  has  the  capa- 
bility to  input  additional  loads  on  elements  and 
nodes . 


EXPLOT  (EXAM  PLOT)  provides  plots  of  the  finite  ele- 
ment structural  model  generated  by  EXAM.  The  plots 
are  isometric  or  two-dimensional  projections.  They 
can  indicate  nodal  points,  elements  and  freedom 
patterns.  Plots  of  any  or  all  of  the  structural 
portions  of  the  vessel  can  be  made. 

LOADER  calculates  the  statically  consistent  nodal 
point  loads  from  the  element  pressures  provided  by 
EXAM  and  from  the  structure's  own  weight.  Addition- 
al dynamic  forces  due  to  ship  motion  are  calculated 
within  the  program.  LOADER  also  calculates  the  out- 
of-balance  forces  for  all  loading  conditions,  and 


the  coefficients  necessary  to  compute  a set  of 
pseudo-inertia  forces  necessary  to  balance  the 
structural  model.  The  out-of-balance  forces  are 
small  in  magnitude,  since  the  vessel  has  been  balan- 
ced in  SHIPMOM,  but  some  small  re-balancing  is  needed 
to  avoid  artificial  reaction  forces  from  developing 
at  the  points  of  rigid  body  support. 

LOADER  also  takes  the  EXAM  output  and  rearranges  it 
in  a manner  suitable  for  input  to  the  DAISY  program. 

7.  CATCH  edits  the  data,  allowing  the  user  to  make  modi- 
fications to  the  LOADER  output  before  input  into  the 
DAISY  program.  The  user  is  able  to  make  changes  to 
any  previously  defined  input  data  by  means  of  addi- 
tional deletions  or  modifications. 

DAISY 

DAISY  (Displacement  Automated  Integrated  SYstem)  is  a 
general  purpose  program  which  performs  linear  elastic  structural 
analyses  of  two-  or  three-dimensional  structures  of  almost  any 
degree  of  complexity  under  statically  applied  generalized 
forces  and  thermal  loads. 

The  DAISY  program  has  a library  of  fourteen  elements 
which  can  be  divided  into  six  categories:  bars,  beams,  membrane 
plates,  bending  plates,  solids  and  substructures. 

Some  special  features  of  DAISY  are: 

1.  Inputs  to  DAISY  can  be  in  standard  fixed  formats  or 
in  any  format  desired  by  the  user.  The  user  may  also 
interface  his  own  preprocessor  and  postprocessor 
programs  with  DAISY. 

2.  Intermediate  results  such  as  element  stiffnesses, 
assembled  stiffness  matrices  and  loads,  can  be  output 
on  magnetic  tape  or  disk  file.  Similarly,  computed 
results  such  as  displacements  and  stresses  can  be 
output  on  magnetic  tape  or  disk  file  for  further 
processing . 

3.  DAISY  automatically  generates  the  total  number  of 
degrees  of  freedom  at  each  nodal  point,  exclusive  of 
rigid  body  supports,  based  on  the  elements  connected 
at  each  node . 

4.  DAISY  is  capable  of  generating  planes  of  symmetry 
and/or  antisymmetry. 
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5.  DAISY  has  a RESTART  capability  designed  to  break 
down  long  computations  into  a series  of  shorter 
runs . 


POSTPROCESSOR  SYSTEM 

The  standard  postprocessor  system  consists  of  DAISYOUT, 
a general  purpose  output  program  which  prints  the  results 
calculated  by  DAISY  with  the  selective  output  of  element 
stresses  and  strains,  and  DAISYPLOT , which  produces  plots 
of  nodal  displacements  and  stress  contours.  The  output  data 
file  produced  by  DAISY  may  print  all  or  part  of  the  results. 
Selective  output  may  be  printed  listing  all  elements  which 
are  stressed  above  certain  user-prescribed  limits.  Normal 
stresses,  shear  stresses  and  stress  concentrations  (Hencky- 
von  Mises  stresses)  are  calculated  and  printed. 


FIGURE  A-l  FLOW  CHART  OF  ABS/DAISY  FINITE  ELEMENT 
STRUCTURAL  ANALYSIS  SYSTEM 


LATERAL  MOMENT  AMPLITUDE,  ft. -ton*  VERTICAL  MOMENT  AMPLITUDE,  ft. -ton* 

WAVE  AMPLITUDE  ft:  WAVE  AMPLITUDE  FC. 


Model  test  data  (From  Reference  Cll}) 


O 25  Kt.  Heavy  Displacement 
\J  30  Kt.  Heavy  Displacement 
A 25  Kt.  Light  Displacement 
& 30  Kt.  Light  Displacement 


Calculated  Results: 


WAVE  LENGTH/SHIP  LENGTH 

FIGURE  A-2  RAO  OF  MIDSHIP  VERTICAL  WAVE  BENDING  MOMENTS  180*  HEADING 


Model  test  data  (From  Reference  Cll}) 


WAVE  LENGTH /SHIP  LENGTH 

FIGURE  A- 3 RAO  OF  MIDSHIP  LATERAL  WAVE  BENDING  MOMENTS  60*  HEADING 
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APPENDIX  B - STRAIN  GAGE  SENSORS  INSTALLED  ON  THE  SL-7 
CONTAINERSHIP 
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TABLE  c - 4 FINE  HESH  15  • PORT 

,Con|il»'  STRESS  : KG/CM 

SELECTIVE  STRESS  OUTPUT  TOR  LOAD  CASE  50  ' 


